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It is shown that a massive Abelian vector boson field can generate the curvature perturbation in 
the Universe, when coupled non-minimally to gravity, through an RA 2 coupling. The vector boson 
acts as a curvaton field imposing the curvature perturbation after the end of inflation, without 
generating a large-scale anisotropy. The parameter space of the model is fully explored, obtaining 
the relevant bounds on the inflation scale and the decay constant of the vector curvaton. 
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Observations provide strong evidence that the Uni- 
verse underwent a phase of inflation in its early history. 
One of the most important consequences of inflation is 
the generation of the curvature perturbation, which is 
necessary for structure formation and observed through 
the CMB anisotropy Quantum fluctuations of suit- 
able fields give rise to a flat supcrhorizon spectrum of 
crturbations through the process of particle production 
Under certain circumstances these perturbations can 
create the curvature perturbation of the Universe. So far 
only scalar fields have been employed for this task. Re- 
cently, however, it has been shown that Abelian gauge 
fields can also work [H, 0|- Indeed, in Ref. Q it was 
shown that, if a vector field obtains a flat supcrhorizon 
spectrum of perturbations during inflation, it can act as 
a curvaton field [|| provided, at some point after infla- 
tion, its mass-square becomes positive and bigger than 
the Hubble scale. In this case the vector field condensate 
oscillates coherently, behaving as pressureless isotropic 
matter Q. Thus, it can dominate the radiation back- 
ground without introducing significant anisotropy, im- 
posing thereby its own curvature perturbation according 
to the curvaton mechanism Hence, the mechanism of 
vector curvaton appears to work using a massive Abelian 
gauge field provided an approximately scale-invariant su- 
pcrhorizon spectrum of its perturbations is created dur- 
ing inflation. In Ref. Q it was shown that this can be 
achieved if the effective mass of the vector field during 
inflation is ~ —2H 2 , where H is the Hubble scale. 
In this letter we investigate a model with an Abelian 
massive gauge field, which is non-minimally coupled to 
gravity, such that the above condition can be satisfied 
during inflation. 

Consider the Lagrangian density: 



the spatially flat FRW metric, which suggests 



R = -6 - + -j = -6(i? + 2H Z ) 



(2) 



where the dot denotes derivative with respect to 
the cosmic time t and H = a/ a, with a being the 
scale factor of the Universe. During (quasi)de Sit- 
ter inflation H~ constant and R ~ — 12H 2 . This 
means that the effective mass of our vector field is 
m 2 ff ~ m 2 — 12aH 2 ~ constant. 

Now, inflation also homogenises the vector field. Fol- 
lowing Ref. [||, we can calculate the spectrum of super- 
horizon perturbations for the vector field. We find that 
the dominant contribution to the power spectrum of the 
vector field perturbations is 
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where k <C aH is the comoving momentum scale and 
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The scale dependence of the power spectrum can be 
parametrised in the usual way as Va k ns ~ , so that 
n s = 1 corresponds to a flat spectrum. Comparing this 
with Eq. (J3j> we find that the spectral index is 



n s — 1 = 3 — 2v =>• n= 



to-4(-j . (5) 



To obtain a scale-invariant spectrum of vector field per- 
turbations we need 



C 



l A„A v + \aRA^A\ (1) 



a^\[l + \{m/Hf] 
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where F, 



fJ,V 



d^A h 



d u A^ is the field strength tensor, 
m is the bare mass of the gauge field and R is the Ricci 
scalar, with a being a real coupling constant. Wc assume 
that a phase of inflation during the early Universe inflates 
away its spatial curvature. In this case we can employ 
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Hence, we see that we need a > ^. If m > H then scale 
invariance is attained only when a is tuned according to 
Eq. However, if m <g; H then scale- invariance simply 
requires a ~ jr. In the latter case m and H do not have 
to balance eachother through the condition in Eq. ^ 
and can be treated as free parameters. We feel that this 
is a more natural setup, so, in the following, we assume 
a « i unless stated otherwise. Since the latest observa- 

6 

tions deviate from exact scale invariance, a should not 
be exactly equal to 1/6. Indeed, according to the 5-year 
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WMAP results n s = 0.960 ± 0.014 at 1-a Q- This im- 
plies that, when m <C H, we need 6a = 1.03 ± 0.01. 

To study the evolution of the vector field we consider 
that, for a homogeneous massive Abelian vector field the 
temporal component A t is zero Q , while the spatial com- 
ponents satisfy the following equation of motion 



A + HA 



±R)A = 0, 



(7) 



where we assume that the homogeneous vector field lies 
along the z-direction with A^ = (0, 0, 0, A(t) ). During 
and after inflation, it is easy to show that 



R = 3(3w - l)H 2 



(8) 



where w is the barotropic parameter of the Universe: 
w ps — 1 [w = |] {w = 0} during (quasi)de Sitter infla- 
tion [radiation domination] {matter domination}. Using 
the above and considering m <C H we can obtain the fol- 
lowing solution for the zero-mode of the vector field 



A = W a + Cci? 



i(3u;-l) 



(9) 



where Wo and C are constants of integration. Thus, the 
growing mode for the vector field, in all cases, scales as 
A (x a. This can be understood as follows. 

As discussed in Refs. [!, H| A^ is the comoving vec- 
tor field; with the Universe expansion factored-out. The 
spatial components of the physical vector field, in a FRW 
geometry are Wi = Ai/a where i = 1,2,3. This can be 
understood just by considering the mass term in Eq. (fTJ), 
which can be written as 



l A^ = WiAt-AAJa 2 ), 



(10) 



where Einstein summation is assumed. Since the La- 
grangian density is a physical quantity we see that 
the spatial components of the physical vector field 
are Wi = Ai/a. Writing the physical vector field as 
Wf, = (0, 0, 0, W(t) ) with W = A/a, we can obtain its 
equation of motion from Eq. ([7]) as 



W + 3HW + m 2 W = , 



which is identical to the one of a massive scalar field and 
we used Eq. (JSJ. When m -C H Eq. (TTT]) has the solution 



W = W +Ca^ w - 1 \ 



(12) 



where Wq and C are constants of integration, consistent 
with Eq. Thus, as long as m <C H, the physical vec- 
tor field develops a condensate which remains constant 
W ~ Wq. This is the physical interpretation of A cx a. 

We can follow the evolution of the vector field conden- 
sate by considering the energy momentum tensor, which 
can be written in the form 



T v = diag(pA, -p±, -Pjl, ~P\\ ), 



where 



PA 



W 2 



±m 2 W 2 



(13) 



(14) 



and the transverse and longitudinal pressures are Q 
P± = UW 2 -m 2 W 2 ) + U2HW+HW + 3H 2 W)W 



(15) 



P\\ 



^ 2 -m 2 W l )-\{2RW+EW+m 2 W)W. 



Thus, the energy-momentum tensor for the homogeneous 
vector field is, in general, anisotropic because p\\ =/= p±. 
This is why the vector field cannot be taken to drive 
inflation, for if it did it would generate a substantial 
large-scale anisotropy, which would be in conflict with 
the isotropy in the CMB. Therefore, we have to investi- 
gate whether, after inflation, there is a period in which 
the vector field becomes isotropic (i.e. p± w py) and can 
imprint its perturbation spectrum onto the Universe. 

Considering the growing mode in Eqs. §§§ and fll2[) . 
from Eqs. (fT4]) and (fl"5]) we see that, during and after 
inflation, when m <C H , we have 



PA 



^m 2 W§ and p A 



\(l-w)H 2 W 2 . (16) 



Hence, the density of the vector field remains roughly 
constant, while the vector field condensate remains 
anisotropic during the hot big bang. 

The above are valid under the condition m <C H ■ How- 
ever, after the end of inflation H{t) oc t , so there will 
be a moment when m ~ H. After this moment, due to 
Eq. ([5]), the curvature coupling becomes negligible and 
the vector field behaves as a massive minimally-coupled 
Abelian vector boson. As shown in Ref. when m > H 
a massive vector field undergoes (quasi)harmonic oscil- 
lations of frequency ~ m, because the friction term in 
Eqs. (7} and (fTTj) becomes negligible. In this case, on 
average over many oscillations, it has been shown that 

W 2 « rr^W 2 0. Hence, Eqs. $3$ and ^ become 



PA — m 2 W 2 and 



\mH [l + |(1 - w){H/m)} W 2 . 



(17) 



(11) The effective barotropic parameters of the vector field are 
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where w± =p±/pa an d Wh =P\\/pa- By virtue of the 
condition m ^ H , we see that, after the onset of the os- 
cillations, u>x,ui|| — > 0. This means that the oscillating 
massive vector field behaves as pressureless isotropic mat- 
ter, which can dominate the Universe without generat- 
ing a large-scale anisotropy. Moreover, its density can be 
shown to decrease as pa oc a~ 3 (like dust) as expected 0]. 
Thus, if the Universe is radiation dominated, pa/ p oc a 
while oscillations occur, so the field has a chance to domi- 
nate the Universe and imprint its curvature perturbation 
according to the curvaton scenario [f|. 
At the onset of the oscillations we have 
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mp 



(19) 
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where we used the flat Friedman equation p = "im 2 p H 2 
with m P = 2.4 x 10 18 GeV being the reduced Planck 
mass. To avoid excessive anisotropy the density of the 
vector field must be subdominant before the onset of os- 
cillations, which means that Wq < mp. 

Let us assume that inflation is driven by some inflaton 
field, which after inflation ends, oscillates around its VEV 
until its decay into a thermal bath of relativistic parti- 
cles at reheating. In this scenario the Universe is matter 
dominated (by inflaton particles) until reheating. Using 
the above findings we can estimate the Hubble scale when 
the vector field dominates the Universe as 



H, 



dom 



mmjm, 1 [ 



(20) 



where T is the decay rate of the inflaton field. If inflation 
gives away directly to a thermal bath of particles then 
we have prompt reheating and T — > where -ff* is the 
Hubble scale of inflation. There is a chance, however, 
that the vector field itself decays before it dominates the 
Universe while still being able to act as curvaton. In this 
case, the density ratio of the vector field at decay is 



n, 



i{ra,r} 
T A 



1/2 



mp 



(21) 



where Fa is the vector field decay rate. 

According to the curvaton scenario the gauge invariant 
comoving curvature perturbation is [f| 



(22) 



where £4 is the curvature perturbation attributed to the 
curvaton field. In a foliage of spacetime of spatially flat 
hypersurfaces Q 



( A = —H — = - — 

Pa 3 pa 



(23) 



dec 



where we used that the vector field decays after the on- 
set of the oscillations in which case pa oc a~ 3 . Note that, 
since (a is determined by the fractional perturbation of 
the field's density, which is a scalar quantity, the pertur- 
bation C,a is scalar and not vector in nature. 

Now, since Eq. (fTTj) is a linear differential equation, 
W and its perturbation SW satisfy the same equation of 
motion. Therefore, they evolve in the same way, which 
means that 8W/W remains constant, before and after 
the onset of oscillations. As shown in Ref. [||, during the 
(quasi)harmonic oscillations of the massive vector field, 
'W 2 , where W is the amplitude of the oscillating 
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physical vector field. From the above we obtain 



2 


SW 




2 


SW 




2 


SW 


3 


w 


dec 


3 


w 


osc 


3 


w 



(24) 



where 'osc' denotes the onset of oscillations and the star 
denotes the time when cosmological scales exit the hori- 
zon during inflation. 



If to <C H during inflation the physical vector field 
(not being conformally invariant) undergoes particle pro- 
duction and obtains an approximately flat superhorizon 
spectrum of perturbations, as shown. Indeed, under the 
condition in Eq. ©, v « f and Eq. © gives @,I3 



aH/2n 



H/2-k, (25) 



i.e. given by the Hawking temperature for de Sitter 
space, exactly as is the case of light scalar fields @j. 
Hence, from Eqs. (|24|) and (|25|) we can write 







3ttW 

Thus, from the above and Eq. (|2"2")) we obtain 



(26) 



(27) 



Using this, Eqs. (|20|) and (|2Tj) . after some algebra, we get 



H*_ ^ C f max{iJ dom ,rA} \ 
mp V^dcc V min{m,r} / 



1/4 



(28) 



The Hot Big Bang has to begin before nucleosyn- 
thesis (which occurs at temperature Tbbn ~ 1 MeV) . 
Hence, max{77dom, Fa} ^ T^ BN /mp. Using this and also 
min{m,r} < H*, we obtain the bound 

H. > C 4/5 ^ 5 (71bn™ 3 p) 1/5 => Vl IA > 10 12 Gc^9) 

where we used that fidec S 1 an d C = 4.8 x 1CU 5 from 
COBE observations. This is similar to the case of a scalar 
field curvaton Q. 

Another bound on the inflation scale is obtained by 
considering that Ta ~ ff 2 w, where g is the vector field 
coupling to its decay products, for which g > m/mp due 
to gravitational decay Thus, maxjiTdom, Fa} > g 2 m. 
Combining with Eq. (|28p we obtain the bound 

K 1/4 > 10 11 GeV , (30) 



H^C^impm) 1 / 2 



where we took m > 1 TeV. 

Finally, an upper bound on inflation scale can be 
obtained by combining Eq. (|27p with the requirement 
Wq < mp, thereby finding 



K 1/4 <10 17 GeV. 



(31) 



where we considered that fldcc ^ 10 2 , in order to avoid 
excessive 11011- Gaussianity in the CMB 0. 

We also need to consider the hazardous possibility of 
the thermal evaporation of the vector field condensate. 
Were this to occur, all memory of the superhorizon spec- 
trum of perturbations would be erased. Considering that 
the scattering rate of the massive vector bosons with the 
thermal bath is r sc ~ g A T we can obtain a bound such 
that the condensate does not evaporate before the vec- 
tor field decays. Since Fgc/r^ oc a -1 , we need to en- 
force this bound at the onset of the oscillations, when 
Tsc ~ j^mpm. Hence, the range for g is 



— <9<( — 
mp \Tnp 



1/4 



(32) 
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where the lower bound is due to gravitational decay. Note 
that, in the case when the vector curvaton dominates the 
Universe before its decay the condensate may not evap- 
orate even if the above upper bound is violated. This 
is because, after domination, the density of the ther- 
mal bath is exponentially smaller than pa by a factor 
of {H dom /H) 2 / 3 . Moreover, even if it does evaporate the 
condensate has already imprinted £a onto the Universe at 
domination rendering the evaporation bound irrelevant. 

The above lower bounds on can be substantially 
relaxed by employing the so-called mass increment mech- 
anism according to which, the vector field obtains its 
bare mass at a phase transition (denoted by 'pt') with 
m/Hpt ^> 1. The mechanism was firstly introduced for 
the scalar curvaton in Ref. 0] and has been already im- 
plemented in the vector curvaton case in Ref. [4[ . 

To illustrate our findings let us consider a spe- 
cific example. Let us choose m ~ 10 TeV and also 
~ 10 -10 GeV such that the temperature at the vec- 
tor field decay is T dcc ~ 10 TeV. Such a particle may 
be potentially observable in the LHC. These values sug- 
gest g ~ 10 -7 , which lies comfortably within the range 
in Eq. (|3"2")) . For the decay rate of the inflaton let 
us chose r ~ 10~ 2 GeV so that the reheating temper- 
ature satisfies the gravitino overproduction constraint 
T Te h ~ \JmpY ~ 10 s GeV. Assume at first that the 
vector curvaton decays before domination Ta > H dom . 
Then Eq. |g8]) reduces to H*/m P ~ lO^C/V^W - Us- 
ing this and Eq. ([57]) we get W /m P ~ 10~ 2 V^dcc- 
Hence, the lowest value for the inflation Hubble scale 
is H* > 10 12 GeV. It can be readily checked that the 
bound in Eq. (|30j) is weaker by a factor 10 -5 . Suppose 
now that the vector curvaton dominates before its de- 
cay T A < H dom . Using Eq. (J20]) we get W > 10~ 2 mp, 
while Eq. 1(25]) suggests ~ (Wo- Taking into account 
the bound Wo < mp, we find that the maximum value 
for the Hubble scale is < 10 14 GeV. The relation be- 
tween and Wo in both cases is depicted in Fig. [TJ 
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FIG. 1: Parameter space for H, and Wo in our example. 

Let us consider now the case when a 96 i. If a = 0(1) 
then, according to Eq. ([6]), a scale invariant spectrum is 
possible only if m~ -H*. Hence, the oscillations begin 
immediately after the end of inflation. With this in mind 
the previous analysis remains valid. In particular, the 
bound in Eq. (|29[) remains the same. However, the bound 
in Eq. (|30|) becomes much more stringent: 

H* > C 2 m P =► V* /4 > 10 14 GeV . (33) 



Hence, in view of Eq. (f3Tj) . we see that the inflation en- 
ergy scale is constrained near that of grand unification. 

In summary, we have discussed a concrete model which 
generates the curvature perturbation in the Universe 
with a single massive Abelian vector boson field non- 
minimally coupled to gravity through an RA 2 coupling. 
The vector field can act as a curvaton, imposing its 
scalar perturbation spectrum well after the end of in- 
flation without introducing a large-scale anisotropy. We 
have shown that there is ample parameter space for the 
model to work by considering all relevant constraints in 
the cosmology. The VEV of the vector curvaton is zero, 
which means that it does not violate Lorentz invariancc 
in the vacuum. Our model does not need to rely on scalar 
fields at all since inflation might take place due to purely 
geometrical effects, such as in /(JJ)-gravity models [l(| 
(e.g. i? 2 -inflation [HI). The remaining challenge is to 
realise our mechanism in the context of a realistic setup 
beyond the standard model [T3 |. 

Recently, vector fields have been employed to drive 
inflation [7| (see also Ref. [IH). To avoid a large-scale 
anisotropy the authors of Ref. [7J] introduce a large num- 
ber of vector fields randomly orientated in space. How- 
ever, they do not consider the generation of curvature 
perturbations, which could proceed along the lines of this 
work, albeit introduced during and not after inflation. 
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